Spin Drag and Spin-Charge Separation in Cold Fermi Gases 
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Low-energy spin and charge excitations of one-dimensional interacting fcrmions are completely 
decoupled and propagate with different velocities. These modes however can decay due to several 
possible mechanisms. In this paper we expose a new facet of spin-charge separation: not only the 
speeds but also the damping rates of spin and charge excitations are different. While the propagation 
of long- wavelength charge excitations is essentially ballistic, spin propagation is intrinsically damped 
and diffusive. We suggest that cold Fermi gases trapped inside a tight atomic waveguide offer the 
opportunity to measure the spin-drag relaxation rate that controls the broadening of a spin packet. 
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Quantum many-body systems of one-dimensional (II?) 
interacting fermions have attracted an enormous interest 
for more than fifty years Contrary to what happens 
in dimensionality D = 2 or D = 3, these systems can- 
not be described by the Landau theory of normal Fermi 
liquids 0]- The appropriate paradigm for ID interacting 
fermions is instead provided by the "Luttinger liquid" 
concept introduced by Haldane in the early 1980s 
The distinctive feature of the Luttinger liquid is that its 
low-energy excitations are collective oscillations of the 
charge or the spin density, as opposed to individual quasi- 
particles that carry both charge and spin. This leads 
immediately to the phenomenon of spin-charge separa- 
tion [H, i.e. the fact that the low-energy spin and charge 
excitations of ID interacting fermions are completely de- 
coupled and propagate with different velocities. 

Recently Recati et al. [1] have proposed to use ID 
two-component cold Fermi gases Q to study spin-charge 
separation. In the case of atoms "spin" refers to two 
internal (hyperfine) atomic states and "charge" to the 
atomic mass density. Kecke et al. Q have pointed out 
that a giant increase of the separation between charge 
and spin modes occurs close to the edge of a harmonic 
potential trap. KoUath et al. "T] have performed a time- 
dependent density-matrix renormalization-group study 
of spin-charge separation in the ID Hubbard model, giv- 
ing quantitative estimates for an experimental observa- 
tion of spin-charge separation in an array of atomic wires. 

In this Letter we expose a new aspect of spin-charge 
separation: namely, spin excitations are intrinsically 
damped at finite temperature, while charge excitations 
are not. The physical reason for this difference is easy to 
grasp. In a traveling spin pulse the up-spin and down- 
spin components of the current are always equal and op- 
positely directed, so that the charge density remain con- 
stant. The relative motion of the two components gives 
rise to a form of friction known, in electronic systems, as 
spin Coulomb draf [1, S, [H El . Of course, no such ef- 
fect is present in the propagation of charge pulses, which 
are therefore essentially free of diffusion and damping, at 



least in the long wavelength limit. By contrast, a density 
pulse in a normal Fermi liquid is always expected to de- 
cay into electron-hole pairs - a process known as Landau 
damping In the following we present a quantitative 
study of this new facet of spin-charge separation, and 
suggest that the theory could be experimentally tested 
in cold Fermi gases. 

We consider a two-component Fermi gas with N atoms 
confined inside a tight atomic waveguide of length L 
(an "atomic quantum wire" ) along the x direction, re- 
alized e.g. using two overlapping standing waves along 
the y and the z axis as in Ref. [5|. The atomic waveg- 
uide provides a tight harmonic confinement in the y — z 
plane characterized by a large trapping frequency uj± ~ 
2tt X 10 kHz [13] ■ The two species of fermionic atoms 
are assumed to have the same mass m and different spin 
cr, (7 =1 or I. The fermions have quadratic dispersion, 
Ek = ?i^fc^/(2m), and interact via a zero-range s-wave 
potential v{x) = gii)S{x) (l5l |. 

The effective ID coupling constant gir, is equal to the 
Fourier transform of the interaction, Vq, and can be tuned 
by using a magnetic field-induced Feshbach resonance be- 
tween the two different spin states to change the 3D scat- 
tering length asD ■ In the limit osd <C a± , where a±_ = 
^f? I {mLu±), one finds gun = 2h? a^Yi / {'ma\) [l^. In the 
thermodynamic limit [N, L oo, N/L — n) the proper- 
ties of the system are determined by the linear density n, 
by the degree of spin polarization C, = (iV| — Ni)/N, and 
by the effective coupling gu). The ground-state energy 
(per atom) £(n, C^id) can be accurately found by solv- 
ing a system of Bethe-^nsate coupled integral equations 
(see e.g. Ref.0)- For future purposes it will be useful to 
introduce the dimensionless interaction 7^ = ravq/ (h^n) 
(for Vq = (71D this quantity coincides with the usual 
dimensionless Yang parameter 7). We also introduce 
the Fermi wave vector fcp — 7m/2, the Fermi velocity 
wf = hk^/m, and the Fermi energy ep = efcp. 

The dynamics of density and spin oscillations is con- 
trolled by the density-density and spin-spin linear re- 
sponse functions, Xppil-,^) and Xs^sA^t^)^ which can 
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be conveniently expressed in terms of the symmetric and 
antisymmetric dynamical local- field factors G±{q,u!) 
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Here xoiQj ^) is the noninteracting ID Lindhard response 
function, which in the long wavelength limit has the form 
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iy{0) = Am/ {■K^'h^n) being the ID density of states. 
The long-wavelength limits of the static local-field fac- 
tors G±{q,Q) are fixed by the compressibility and spin- 
susceptibility sum rules, i.e. 



(3) 



where n = n^'^{d'^[ne{n,Q,giYy)]/din?}^^ and Xa = 
n[9^£(rt, C, 5'id)/9C^Ic=o]~^ are the compressibility and 
the spin susceptibility of the interacting system, kq = 
v{Q)/n?' and xm — ^{^) a^'^ the same quantities for the 
noninteracting system. 

At finite frequency, however, the real and imaginary 
parts of G-{q,uj) diverge, in the long-wavelength limit, 
as io^ /{q^Vq) and uj/{q'^Vq) respectively [3]. No such di- 
vergence exists in the density channel. As a result, the 
small-g behaviors of the density-density and spin-spin re- 
sponse functions are dramatically different. Indeed, using 
Eqs. (HI)-© and the known form of the singularity in G- 
it is easy to show that the inverses of these functions have 
the following form: 
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where, following the nomenclature introduced in Ref. [lOj, 
we have introduced the "spin mass" m„ and the inverse 
of the spin-drag relaxation time - both functions of 
uj and temperature T (23 |. 

We now analyze in detail the physical implications of 
Eqs. (gl) and p. 

According to Eq. ([4]) the density-density response func- 
tion has an acoustic pole at tu = Vpq where Vp = v'^kq/k. 
This result is easily recognized to be in agreement with 
the bosonization result Vp = vpKpKo/ k due to the 
relation Kp = v-p/vp that holds for a Galileian invariant 
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FIG. 1: (Color online) Spin-drag relaxation rate (in units 
of ef/A and divided by 7'^) as a function of the reduced tem- 
perature kBT/{2eF). In the inset we show a zoom of the 
low-temperature region < ksT /{2eF) < 1: the filled circles 
represent the analytical result ([8]). 



system. The speed of sound Vp has the following behav- 
iors, Vp = i'f(1 + 7/71"^ + ...) in the weak coupling 7 — > 
limit, and Vp = 2wf(1 — 41n2/7-f...)in the strong cou- 
pling 7 +00 limit 01 ■ Eq. ^ docs not incorporate 
the damping of this acoustic mode associated with the 
quadratic energy dispersion of the particles: in fact, this 
damping appears only at higher orders in q and goes to 
zero as q^ ITI ]. 

According to Eq. ([5]), the spin-spin response func- 
tion has an acoustic pole at a; = v„q where = 
v'^{'m/m„)Xcro/X(y The bosonization result for the spin 
velocity is v„ = vp K a-Xa-o / Xa- Hi with Kg. — 1 due to 
spin-rotational invariance. The two results for Wg. coin- 
cide only if the following non-perturbative relation be- 
tween the spin mass and the spin susceptibility holds: 
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In the weak coupling limit = vp{l — 7/7r^ + ...). In 
the strong coupling limit the spin velocity goes to zero as 
Va = 2'k'^vf/(3j) + ... [1]. From these limiting behaviors 
we find that mfj/m = 1 -t- 7/7r^ -I- ... in the weak coupling 
limit and that the spin mass diverges linearly at strong 
coupling, TOcr/m = 37/(27r^) -I- .... 

Let us now examine the spin-drag relaxation time, 
which is responsible for the damping of the spin mode. In 
the unpolarized case iV-f = and within second-order 
perturbation theory the spin-drag relaxation rate (at zero 
frequency) is given by the formula [§] 
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FIG. 2; (Color online) Space-time evolution of a Gaussian 
spin packet of initial width w = lOfcp^ in the ballistic limit. 
s(x,t) (in units of fcp^) is shown as a function of x (in units 
of kp^) and t [in units of h/{2£F)]. In this example we have 
chosen 7 = 0.6. The left and right components of the packet 
propagate at essentially the Fermi velocity and their width 
does not change. 



Here xo{q,^,T) is the ID Lindhard function at finite 
temperature In Fig. [T] we plot t~^{T)/^'^ as a func- 
tion of temperature. We clearly see that the s pin -drag 
relaxation rate goes to zero linearly for T — > [2l|. In- 
deed, using Eq. ([7]) it is possible to show that 
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With increasing temperature the inverse spin-drag relax- 
ation time first saturates and then decays to zero rather 

slowly: t;^{T) IQ^-^^-f'^lkBT / {2eF)]-^/^eF /h, for 

Vq = 5iD- 

As discussed in Refs.|jand|7|, localized spin and density 
packets can be created by short off-resonant and state 
selective laser light pulses focused to a spot size w such 
that R ^ w ^ k^^ , R being the size of the atom cloud 
and fcp ^ the average interatomic distance. 

The time evolution of these packets is deter- 
mined by the equations ^^^-'^(g, cj)n((7, cj) = and 
Xg^^g^(g, w)s(<7,u;) = 0. Making use of Eqs. (g]) and ^ 
these can be readily converted into two second-order par- 
tial differential equations for the density and the spin 
density: 
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FIG. 3: (Color online) Space-time evolution of a Gaussian 
spin packet of initial width w = lOfcp ^ in the presence of a 
spin-drag relaxation time Tad ~ 0.04/i/eF (corresponding to 
7 = 0.6 and T = £f/A;b). Notice the spreading and attenua- 
tion of the packet as time progresses. 



this equation with initial conditions n{x,0) = no(x) and 
dtn{x,t)\t=o = is given by the d'Alembert formula 
n{x,t) = [nQ{x + Vpt)-\'nQ{x — Vpt)]/2. The time evolution 
of a density packet is therefore quite simple: the packet 
splits into a left-moving and a right-moving component, 
each one preserving the shape of the initial profile while 
moving ballistically. 

The second equation is a damped wave equation that 
depends on two parameters, i.e. the spin velocity 
and the spin-diffusion constant D^. In the undamped 
Da — > -foo limit its solution is given by the d'Alembert 
expression s(x,t) — [so{x + v^t) + sq(x — Vai)]/2. In 
this limit both density and spin packets move ballisti- 
cally (with different velocities). This situation is shown 
in Fig. [2] for an initial Gaussian spin packet sq{x) = 



[(27r) 



-1/2 



W 



^] exp [-x^ /{2w^)] of width w = lOfcp 



For a finite value of the spin diffusion constant the 
dynamics of the spin packet becomes noticeably differ- 
ent from that of a density packet. The solution of the 
damped wave equation satisfying the initial conditions 
s{x,Q) = So (a;) and dts{x,t)\t=() = can be written as 
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{v-^df-d^,)n{x,t) = 
{v-^d^^-dl)s{x,t) + D-^dts{x,t) = Q 



(9) 



where we have introduced the spin diffusion length = 
VDaTsd = Da/va, ujq = ^(2gi^)2 - l/(2Tsd) (the com- 
plex square root is defined here with a positive imagi- 
where we have introduced the spin diffusion constant nary part), and so (g) = dxso{x) exp {—iqx) (Fourier 



Da = vlTsd{T) = nTsd{T) / (rriaXa) (Einstein relation). 

The first equation is a ID wave equation that de- 
pends only on the density velocity Vp: the solution of 



transform of the initial spin density profile). 

In Fig. [3] we show the time evolution obtained from 
Eq. (fTUl) of the same Gaussian spin packet as above for 
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FIG. 4: (Color online) Space-time evolution of a Gaussian 
spin packet of initial width w = 20fcJ^ in the presence of a 
spin-drag relaxation time Tad ~ 0.12fi/£F (corresponding to 
7 = 1 and T = 2eF /fee ) • The packet spreads out by diffusion 
and does not split into left- and right-moving components. 



7 = 0.6 and T — e-p/k-Q. The spreading and attenua- 
tion of the packet are quite noticeable. If the spin-drag 
relaxation time is sufficiently short {vaTsd < w) it may 
even prevent the splitting of the initial profile into two 
peaks: in that case the evolution of the packet is hardly 
distinguishable from ordinary spin diffusion in higher di- 
mension (see Fig. HI . 

In summary, we have shown a new aspect of spin- 
charge separation in one-dimensional Fermi systems. Not 
only spin and charge propagate independently at differ- 
ent speeds ~ they are also damped at different rates. The 
propagation of charge excitations is essentially ballistic - 
any residual damping vanishing in the limit of a smooth 
density profile. On the other hand, spin density exci- 
tations are intrinsically damped, no matter how smooth 
the profile, and undergo diffusion, which may completely 
suppress the ballistic behavior. We believe that two- 
component ID trapped Fermi gases, in which spin pulses 
can be created and monitored at different times, are ide- 
ally suited for an experimental verification of these ideas. 

This work was supported by NSF Grant No. DMR- 
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